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Doping of substrates at desired locations is a key technology for spin-based quantum memory 
devices. Focused ion beam implantation is well-suited for this task due to its high spacial resolution. 

In this work, we investigate ion-beam implanted erbium ensembles in Yttrium Orthosilicate crystals 
by means of confocal photoluminescence spectroscopy. The sample temperature and the post¬ 
implantation annealing step strongly reverberate in the properties of the implanted ions. We find 
that hot implantation leads to a higher activation rate of the ions. At high enough fluences, the 
relation between the fluence and final concentration of ions becomes non-linear. Two models are 
developed explaining the observed behaviour. 
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I. INTRODUCTION 

The development of systems and devices suitable for 
quantum communication protocols has largely advanced 
during the last two decades. These systems are optimised 
for a single application such as data processing, data trans¬ 
mission or data storage. In order to develop a feasible 
quantum-information infrastructure, it is necessary to com¬ 
bine the best properties of those indi vidu al quantum sys¬ 
tems within a so-called hybrid systerrP^. An example of 
a promising hybrid system are nuclear and electron ^ins 
coupled to a superconducting resonators and qubits^nS. A 
great leap towards the realization of such an interface had 
been made in cir cuit QED experiments with a variety of 
spin-ensemble^^I^®^. Spin-ensembles can be prepared lo¬ 
cally on the substrates, thus allowing an arrangement of 
several bit-units with customised properties on one chip. 
This enables a simple and available technological process 
for fabrication of quantum circuits. 

We would like to emphasise the ensembles of rare-earth 
ions (RE) due to their unique properties. The optically ac¬ 
tive 4f-electrons of the RE elements are semi-shielded from 
the external fields by the 5d and 6s orbitals, which results in 
long optical and microwave coherence times. Kramers RE 
ions possess large magnetic moments, which can be used 
for the manipulation of their spin states with microwaves. 
The realization of such spin ensembles by focused ion beam 
(FIB) implantation technique has many advantages, such 
as freedom in positioning of different ensembles in one crys¬ 
tal, variability of their concentrations and design of their 
arrangement in a single mask-less procesJ^. 

One of the crucial requirements for the host substrates 
is a low local symmetry of the surrounding oxygen ions. 


which stabilises the RE ions in the lattic#^*^^. This prop¬ 
erty is inherent in Y2Si05, Y3AI5O12, YAIO3 and YVO4 
crystals. In this work, we present in-depth analysis of the 
FIB-implanted erbium ensembles in Y2Si05 (YSO) crys¬ 
tals. This crystal is known to have the longest measured 
optical coherence time of about 4 ms i nside the telecom C- 
band at a wavelength of about 1.5 However, for¬ 

mation of defects and dislocations during the implantation 
process aggravates the optical properties of the ions. There¬ 
fore, a deep understanding of ion-matter interactions and 
activation rates in a particular substrate is indispensable. 
The first results on the activation rate of the implanted 
Erbium ensembles have already been reported irfi^. 

In this work, we demonstrate a more fundamental anal¬ 
ysis of the focused ion beam implanted Erbium ensembles 
in the Y2Si05 (YSO) crystals in dependence on the prepa¬ 
ration process. The implantation of YSO crystals is per¬ 
formed in a wide range of fluences at different tempera¬ 
tures of the substrate. Optical properties of the implanted 
ions are analysed in dependence on fluence, ion concentra¬ 
tion and ion activation by means of confocal photolumi¬ 
nescence. We compare the results of the photolumisence 
measurements to the results of microwave spectroscopy of 
the implanted sample^^^. 


II. EXPERIMENTAL TECHNIQUES 

Focused ion beam implantation 

Samples are implanted in a focused ion beam (FIB) ma¬ 
chine EIKO-IOO with an accelerating voltage of 100 kV un¬ 
der high-vacuum conditions of 10“^ mBar. Schematic of 
the implantation process is demonstrated in Fig.[^(a). In 
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Figure 1. (a) Schematic of the implantation process. The orientation of the optical symmetry axes is shown accordingly to the 
sample fabrication and direction of the ion-beam. The ion channelling direction coincides with the D2-axis and is perpendicular to 
the implantation direction, (b) Mass-spectrum of the Er"'"^ peak. The distribution of Erbium isotopes is given by bars according 
to their natural abundance. The arrow points to the selected position for the implantation, (c) SRIM-simulated data for the 
ion-distribution (blue) and the vacancy-distribution (red) and RBS measured ion-distribution (green). The RBS result corresponds 
to the fluence of 2 x 10^® cm“^, below which the concentration is under the detection limit. Depth of the ion-maximum is 50 nm for 
both SRIM and RBS results, (d) Scanned luminescence pictures of the implanted patterns, intensity in Counts. Letters a, b, c, and 
d in the middle of patterns correspond to the implantation process as follows: a - implanted at 600 K and annealed; b - implanted 
at 300 K and annealed; c - implanted at 600 K and not annealed; d - implanted at 300 K and not annealed, (e) Level-diagram of an 
Er®+ ion. Excitation is schematically shown with the red arrows. 


the current work, the FIB operates with a Liquid Metal 
Alloy Ion Source (LMAIS) Au 78 . 4 Siii. 6 Erio developed by 
Melnikov et al.l^^. Ions are separated by a built-in Wien 
mass-filter with a resolution of 2 a.m.u. For the Er++ peak, 
Fig.[T](b), it is possible to make a selection of the isotopic 
mixture with the lowest content of ^®^Er (to avoid the spec¬ 
tral broadening from the hyperfine structure). The distri¬ 
bution of the Erbium isotopes relative to the mass spectrum 
peak is shown in Fig.[^(b). In this work, the 90%-pure 
i70Er++ isotope is implanted with a Wien-filter setting as 
shown with an arrow in Fig.[^(b). 

For YSO crystals, the sputter rate is low due to the 
high binding energies of the ion j^^ l ^° [ The distributions 
of ions and induced vacancies are simulated with the SRIM 
softwar^^, Fig.[^(c). As the channelling direction coincides 
with the D2-axis and is perpendicular to the direction of 
the ion-beam, Fig. [i](a), the SRIM simulation is applicable 
for the ion-distribution. From the ion- and defect-profiles, 


we estimate final ion-distribution to be 50 nm deep for the 
ion maximum and 32 nm of the ion straggle. For the sam¬ 
ple implanted with a fluence of 2x10^^ cm“^, a Rutherford 
Back Scattering (RBS) profile is measured and compared 
to the simulation, Fig.[^(c); its maximum is found at 50 nm 
from the surface with the straggle of 70 nm. Thus at higher 
fluences, migration of the ions takes place during the an¬ 
nealing and, therefore, the distribution of the ions broad¬ 
ens. For lower fluences, the RBS is not measured as the ion 
density is below the detection limit. 

To reduce the defect density, we apply two methods: in¬ 
crease of the substrate temperature during the implanta¬ 
tion and variation of the post-implantation annealing pa¬ 
rameters. To increase the substrate temperature during 
the implantation, a special sample-stage is designed with 
a heating element and a temperature detector with a tem¬ 
perature limit of 600 K, which is close to the Debye tem¬ 
perature of Y 2 Si 05 , Td = 580 KP^. Annealing procedures 
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are varied from the rapid thermal annealing (RTA) in N 2 
atmosphere to the long thermal annealing in Argon and 
air. Annealing for several hours in Argon atmosphere or in 
air is found to be the most effective. An example of four 
implanted patterns is given in Fig.[^(d). Other samples are 
implanted with larger patterns of 1 mmx 1 mm in size. 

Measurement technique 

Samples are characterized with 656 nm CW laser to 
achieve the up-conversion excitation at room temperature 
and at 4 K in a confocal setufP^. Erbium luminescence is 
measured in reflection and directed to a grating spectrom¬ 
eter. The detection range is from 490 nm to 650 nm, see 
Fig. UK e). For 4K measurement, sample is placed on a cold 
finger in a cryostat with access to the sample through a 
window. To measure relaxation times, samples are excited 
with pulses of 50 /rs and the following decrease of the in¬ 
tensity is recorded, [C] and Fig.[^ 


III. RESONANCE ENERGY TRANSFER 

At low concentrations ions, the intensity of the lumi¬ 
nescence is linearly proportional to the number of the ions. 
When the concentration increases, the non-radiative energy 
exchange between the ions becomes dominant, and the to¬ 
tal emitted energy distr ibute s between the radiative and the 
non-radiative processe^^lllll Qjjg jg described by 

non-radiative resonant energy transfer (RET), which hap¬ 
pens either due to a perfect match between the absorption 
and the emission energies of the neighbouring ions or with 
assistance of the ion-phonon interactions, which compen¬ 
sate mismatches of the ions energies. This effect is also 
known as a saturation and quenching of the luminescence 
intensity with the increase of concentratiorP^ and is de¬ 
scribed by the Forster lavP^. 

From the RET theory and Forster law, the non-radiative 
energy transfer rate, WjsfR, is due to the electric-dipole- 
electric-dipole (ED-ED) interactions and is proportional to 
R-® 

W^n = C^, ( 1 ) 

where R is the distance between the donor and accep¬ 
tor ions, C is the concentration of the acceptors, Rp is 
the Forster radius, at which energy on the radiative and 
non-radiative exchange is distributed equally. As we char¬ 
acterize optical transitions of Erbium, which are electric- 
dipole allowecP^l^^, we consider energy transfer only due 
to the ED-ED interaction and neglect the higher order in¬ 
teractions and the interactions involving magnetic dipoles 
(MD). To apply the RET to our experimental data, we 
need to integrate the energy transfer probability for one 
donor-acceptor pair over the possible neighbours in the 
volume2IH2Sl: 

F = C'y J j ^dr^drydr^, (2) 

with R^ = The ion-concentration varies over 

the implanted layer, Fig.[^(c), thus the interatomic dis¬ 


tance R varies as well. To simplify, we introduce a mean 
interatomic distance, Rm, which is derived from the im¬ 
planted fluence F and the ion straggle d: Rm = ^/d/F, 
and an average concentration Cav = F/d, which is the 
concentration at half-maximum of the ion-distribution, see 
Fig. UK c). This average concentration is twice smaller than 
concentration at the maximum of the ion-distribution peak. 
Thus, we obtain the following expression for the RET rate: 

1 47r 

F = Ca. R% ^d^r = (3) 


The fit-function for the luminescence intensity on the flu¬ 
ence in this case has the following form: 
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where S is area of the measured spot and /q is intensity per 
ion. 

At low ion-concentrations, RET is absent or negligibly 
small, and the luminescence intensity is proportional to the 
number of ions, and therefore to the fluence: / = Iq S F. 
When the ion density increases, the impact of RET on the 
luminescence yield becomes significant, and a part of the 
energy goes into the non-radiative energy exchange between 
the ions. In the result, the luminescence intensity per ion 
decreases with the number of ions as a function of their 
concentration: Iq = /o/(l + r). This equations are further 
used to explain the experimental dependence of lumines¬ 
cence on fluence. 


IV. EXPERIMENTAL RESULTS 

Shapes of the spectra 

Luminescence study of the annealed and non-annealed 
samples is presented in Fig.[^ The spectral shape of the 
Erbium optical transitions contains information about the 
surrounding crystal field and, consequently, about the ion 
activation in identical lattice sites. This information is 
given by the positions and linewidths of the spectral lines 
and overall shape of the spectra. Variety of local crystal 
fields from site to site results in a broad spectrum. This can 
be seen for the non-annealed samples in Fig. [^(a): For all 
non-annealed samples with the fluences above 10^^cm“^, 
spectra take a broad shape and specific sub-level transi¬ 
tions are not distinguishable. The implantation tempera¬ 
ture is found essential at low fluences, when noticeable ion- 
activation without annealing is favoured. This effect can 
be seen in the normalized spectra in Fig.[^(a) for the flu¬ 
ences below 10^^cm“^ as more distinguished shapes of the 
emission lines of the non-annealed sample. For the room- 
temperature implanted samples, spectra are similarly broad 
at all fluences. 

Annealing time and atmosphere have crucial effect on 
the stabilization of the ions in the lattice positions. We 
found that long annealing in air leads to an optical re¬ 
sponse, which is the most close to the doped as-grown sam¬ 
ple. Comparison of the annealing procedures can be found 
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Figure 2. (a) Luminescence spectra of the '* 83/2 ~^^Ii 5/2 manifold for annealed and non-annealed samples implanted at temperatures 
of T = 300 K and 600 K temperatures, indicated in each graph. All spectra are normalized to their maximal value. For the annealed 
samples, spectra are similar for all flnences. (b-c) Deconvolution of "^ 83/2 —>■ “^115/2 manifold spectra of doped as-grown and implanted 
samples, which are measured at (b) 300 K and (c) 4K. For the spectrum at 4K (c), the independent lines are assigned to the 
transitions to specific sub-levels of the ground multiplet. 

in 1^ Further, we discuss the results for the samples an¬ 
nealed in air. 

For the air-annealed samples, the spectral shapes are in¬ 
dependent of the fluence^S and are similar for both im¬ 
plantation temperatures, Fig.[^(a). At the fluences above 
10 ^®cm“^, several new emission lines appear while a few 
other are enhanced in intensity. This occurs as a result 
of the high density of ions and probable formation of clus¬ 
ters. Spectra of the implanted samples reveal differences 
from the doped as-grown crystal, which are similar over all 
fluences and implantation temperatures. This proves the 
reproducibility and stability of the implantation-annealing 
processes, which nevertheless lead to minor differences in 
local crystal symmetry as compared to the doped as-grown 
crystal. 

Spectral lines 

More detailed analysis of the implanted ions is made for 
the manifold ^ 83/2 —>-^ 115/21 as shown in Fig.l^(b,c). At 
room temperature, positions of the deconvolutea peaks are 
similar for all samples and are most closely described by a 
Lorentzian function. The minimal linewidth of the spec¬ 
tral lines at room temperature is found to be 1.32 meV 
(319.0 GHz) for the implanted samples and 1.29 meV 
(311.7 GHz) for the doped as-grown. The maximum 
linewidth at room temperature is 4.38 meV (1054.5 GHz) 
for the both implanted and doped as grown samples. 

At low temperatures (Fig.[^(c)), only eight transition 
lines remain. The lines labelled 1 and 2 exhibit purely 
inhomogeneous line shape and linewidths of 0.56 meV 
(134.8 GHz) for the implanted samples and of 0.83 meV 
(201.3 GHz) for the doped as-grown. However at longer 


wavelengths, lines become homogeneously broadened, and 
the linewidth increases to 1.53 meV (369.7 GHz) for the im¬ 
planted samples and 1.64 meV (396.3 GHz) for the doped 
as-grown, as can be seen for the lines 5-8 in Fig.[^ All eight 
transition lines are associated with the transitions from the 
lowest level of the "^ 83/2 multiplet to the eight sublevels of 
the ground multiplet "^Ii 5/2 as shown in Fig.^(c). They cor¬ 
respond to the Site 2 of Erbium in YSO noted by Bdttger 
et al.^. 

The variation of the linewidth and the line-shape with 
the wavelength at low temperature is demonstrated in 
Fig.|( a). The very first transition line has a purely in¬ 
homogeneous, Gaussian shape, while for the other lines, a 
homogeneous, Lorentzian, component becomes more rele¬ 
vant and the linewidth increases. 

In order to separate the homogeneous and inhomoge¬ 
neous broadening components and to see how their values 
vary with the wavelength, the low-temperature transition 
lines are fit with a Voigt fit-function. Fig. [^(c) . The derived 
components are plotted in Fig.|^(b). The inhomogeneous 
component remains similar for all eight lines (as demon¬ 
strated with the solid lines in Fig.[^(b)). The homogeneous 
component appears to be significant in the third transition 
line and increases until it is completely dominant in the last 
three lines. 

Appearance of a large homogeneous broadening is ex¬ 
plained with the transitions within the ground-state mul¬ 
tiplet, as shown in Fig.[^(c). After a radiative transition 
to one of the sub-levels of the ground ^Ii 5/2 multiplet, an 
electron relaxes further to the lowest sub-level. Each tran¬ 
sition line is assigned to a particular transition from the 
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"^ 83/2 multiplet to the ground state sub-level as depicted in 
Fig.|( c), and the homogeneous broadening part is found 
to be proportional to the energy separation between the 
higher-energy sub-levels of the ground multiplet and its 
lowest sub-level. 

For the direct transition from “^ 83/2 to the lowest sub- 
level of the ground multiplet at low temperature, when 
phonon modes are “frozen”, the natural linewidth (40 kHz 
or 0.16x 10“® meV) is dominated by an inhomogeneous 
broadening due to the variations in the crystal-field strains 
experienced by each ion, differences in the transition en¬ 
ergies for different isotopes and/or impact from the hy- 
perfine structure. Thus, in case of a larger inhomogeneous 
linewidth, ions are experiencing more differences in the sur¬ 
rounding crystal fields. The first and the second transition 
lines in Fig. K c) have purely Gaussian shape. These lines 
are assigned to the transitions between the lowest state of 
the ^ 83/2 manifold and two lowest states of the ‘^Ii 5/2 man¬ 
ifold, Fig.[^(c). 

As mentioned earlier, the implanted samples show nar¬ 
rower linewidth as the doped as grown sample: 0.56meV 
versus 0.83meV, respectively. The width of the spectral line 
has contribution from hyperfine structure, isotope shifts 
and crystal field. Implanted Erbium has « 90% of ^^°Er- 
isotope which can be the reason for a smaller linewidth. 
However, the inhomogeneous broadening of a transition 
line due to t he pr esence of different isotopes is taken from 
Haynes et al.l^SEZl and equals to 0.02 meV. The broadening 
due to the hyperfine structure of the Er^®^-isotope is ex¬ 
pected to be in the same order as the estimated isotopic 
broadening. Therefore contributions from the hyperfine 
structure and isotopic broadening are neglected. The dif¬ 
ference between the emission linewidths of the implanted 
and doped as-grown samples can be explained by formation 
of a more homogeneous local crystal field. 

The destruction of the substrate by the implantation and 
recovery of it by the annealing creates a more homogeneous 
crystal structure and thus results in a smaller number of 
differences in the crystal-fields. This more homogeneous 
crystal structure is not due to lower number of defects in the 
lattice, but due to a more homogeneous defect distribution. 

However at room temperature, the linewidth in the 
doped as-grown sample is smaller as compared to the im¬ 
planted one. The dependence of the linewidth on the tem¬ 
perature is commonly described by the McCumber- 8 turge 
relatioiP^ 


r(T) = Fo + a 



X® X 

(e“ - !)■ 


-dx. 


( 5 ) 


where Tp is the Debye temperature of the material, Fq is 
the linewidth at 0 K temperature, T is the actual temper¬ 
ature and a is the coefficient proportional to the strength 
of the electron-phonon interactiorPHHini. Thus, the influ¬ 
ence of the phonons on the spectral lines appears in the 
value of a. From the known linewidth values, we derive 
values of a for our samples. For the doped as-grown sam¬ 
ple, a = 93cm“^, whereas for the implanted sample it is 
equal to 155 cm“^. This difference in the coefficients indi¬ 
cates that the ion-phonon coupling is much stronger in the 


implanted samples. The stronger ion-phonon coupling can 
be interpreted as a confirmation of a specific homogeneous 
defect-entourage, which leads to the narrower spectral lines. 

Implantation yield 

Not every implanted ion can settle in a desired lattice 
site. Ratio between the implanted fluence and the final 
amount of the activated ions is described by the implanta¬ 
tion yield. It can be derived from the number of lumines¬ 
cent ions normalized to the total number of the implanted 
ions. We estimate the yield by comparing the luminescence 
of the doped as-grown sample to the implanted sample and 
by taking the ratio of the concentrations corresponding to 
equal intensities, as described irfi^. The measured volume 
is assumed to be equal to 300 nmx 300 nm x 1 /rm, and thus 
the implantation yield is found to be 20 % for the room- 
temperature implantation and 22 % for the hot implanta¬ 
tion. 

Energy transfer between the implanted ions 

For the previous calculation of the implantation yield, we 
plotted the integrated luminescence intensity as a function 
of fluence, Fig.[^ The intensity is integrated over three 
manifolds separately: ‘^F 7/2 —>-"^Ii 5 / 2 , ^Hii /2 —t"^Ii 5/2 and 
'^ 83/2 —>'"‘li 5 / 2 - For each manifold, the dependence is anal¬ 
ysed according to Eq. Q. 

For the annealed samples, two saturation steps are ob¬ 
served, as shown in Fig. |(a ,b). the first saturation step 
appears at a fluence around 10 ^"^ cm“^, after which the in¬ 
tensity continues to increase with the fluence until it is 
completely saturated at fluences of 10^®cm“^. To explain 
these two steps, we consider the limitation of the Erbium 
concentration in YSO for growth, which is approximately 
10% of the total number of Yttrium iona^S and equals 
Cth = 2.4 X 10^°cm“^. This concentration at the max¬ 
imum of the ion-distribution corresponds to a fluence of 
3.84 X 10^^cm“^. If the 10%-limit is exceeded, creation of 
Er-8i-0 formations is observecP^. 

With respect to this, we explain the two saturation steps 
as following: At fluences of lO^'^cm”^, the local volume 
concentration reaches the growth limit of Cth- Ions are 
embedded into the lattice so close to each other that RET 
takes place, and a saturation step is observed. However 
during the annealing, the crystal structure tends to stay as 
Y 28 i 05 crystal without cracks and clusters. This results in 
migration of the implanted ions to the closest vacant posi¬ 
tions created by implantation, which are dislocations in the 
damaged layer (distribution shown in Fig.[^(c)). As a re¬ 
sult, the highest ion concentration remains relatively same 
while the ion-straggle d increases, and therefore the to¬ 
tal luminescence intensity increases with luminescent yield 
per ion (and concentration) remaining relatively constant, 
Fig.|^(f). The diffusion of the ions is limited by the thick¬ 
ness of the available damaged layer, which equals to the 
sum of the depth of ion-maximum, h^ax = 50 nm, and half 
of the ion-straggle, d/2 = 16 nm, see Fig. [^(c). The result 
of Ratherford back scattering for a sample implanted with 
a fluence of 2 x 10^^ cm“^ is presented in Fig. ([b(c)) and is 
the experimental proof of the diffusive behaviour. 

When the avaliable layer appears to be completely filled 
with threshold concentration of Cth, the reached straggle 
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Figure 3. (a) Linewidth values for the eight transition lines of ^ 83/2 ^^Ii 5/2 manifold observed at 4K. thedotted line describes in¬ 
crease of the linewidth value with the wavelength, (b) Lorentzian and Gaussian components of the spectral lines of the ^ 83/2 —>'^Il5/2 
manifold. The Gaussian component remains relatively constant (solid line), while the Lorentzian component rises with the wave¬ 
length (dashed line). The increase of the Lorentzian component and thus the broadening of the spectral lines is suggested due to 
electron relaxation from the higher sublevels to the ground-state sublevel of the ground multiplet '^Ii 5 / 2 , (c). 


dmax becomes constant, and the concentration begins to in¬ 
crease again. With this rise of concentration, luminescence 
quenching is observed again, yielding into a continuation of 
the previous slope. We name this behaviour as Model 1, 
which is schematically demonstrated in Fig.[^(f). 

Experimental data for the annealed samples are fitted 
with a function based on Eq. which is adapted for 
the Model 1: for the fluences below the threshold fluence 
Fth ~ 10^‘^cm“^, the concentration changes proportion¬ 
ally with the fluence; above the Fth-, the concentration 
remains constant until the volume limit is reached, after 
which the concentration increases proportionally with the 
fluence again. Details on the fit functions are given in [D| 
However, above the threshold fluence Fth, the concentra¬ 
tion may still be slowly increasing, simultaneously with 
the ion-straggle, Model 2. In this case, the intensity de¬ 
pends on both the straggle and the fluence, and its slope 
changes. This appears to be the case for the room tem¬ 
perature implantation, Eig.|^(b), while hot implantation is 
well-described with Model 1, Fig. |4](a). 

The solid line in Fig.l^l^b) shows the fit until fluences of 
2 X 10^"^ cm“^ using Eq. Q. By fitting Model 1 to the full 
range of fluences, we obtain good agreement with the data, 
as shown by the dotted line. For the both fits, the Forster 
distance is derived from the fit to be equal to 2.31 nm and 
2.18 nm for the hot and room-temperature implantations, 
respectively. These values are similar to the Forster dis¬ 
tances numerically evaluated in the literature for the tran¬ 
sitions between ^Ii 3/2 and ^Ii 5 / 2 , which are approximately 
1.73 nniil^. The Forster distances are expected to vary 
for different manifolds, however, they should stay in the 
same order of values. The threshold fluence for the hot 
implantation is equal to 1.5 x 10^^ cm“^, and for the room- 
temperature implantations it is equal to 3 x 10^^cm“^. 
Earlier luminescence saturation for the hot implantation 
is related to the higher activation rate of the ions at the 


crystallographic sites during the implantation and, thus, 
higher final concentration at equal fluences. 

Eor the non-annealed samples, we observe only one satu¬ 
ration step at fluences of 10^® cm“^. The energy transfer in 
the annealed samples can be assisted by phonons. There¬ 
fore, it can occur for larger mismatches between the transi¬ 
tion energies. In non-annealed areas, a large variation of the 
crystal-field strengths reduces matching between the emis¬ 
sion and absorption energies of the neighbouring ions. Ad¬ 
ditionally, there are very few phonon energy-states present 
in the damaged crystal lattice, and, therefore, it requires a 
much higher fluence to achieve the luminescence saturation 
effect. 

Eor the non-annealed samples, we obtain the Forster 
distance of 0.54 nm for the room temperature implan¬ 
tation, Fig.|^(d), and 0.61 nm for the hot implantation. 
Fig .[^(c), with the ion straggle d taken from SRIM as 
32 nm, Fig.[^(c). These values are already close to the 
smallest possible distances between the two ions in the lat¬ 
tice, which is approximately 0.5 nm. At such distances, 
luminescence quenching may happen via an electron trans¬ 
fer (Dexter energy transfer)p21. The luminescence intensity 
per ion in the non-annealed samples is higher in the case of 
hot implantation, which is again due to the partial direct 
activation of the ions in the lattice. 

Model 1 and Model 2 are also supported by the dynamics 
of the relaxation time value. The relaxation time shortens 
as the concentration of ions and the number of interactions 
between them increas^^. At the lowest fluences, our mea¬ 
sured relaxation times are equal to 4.2/is and reduce to 
4.0/iS at the fluences of 10^^cm“^. Their values remain 
equal to 4.0 ^s in entire diffusion range. Around the final 
saturation step, the relaxation times decrease and reach 
3.2 fj,s. 

The effects of RET and the concentration non-linearity 
are also observed in circuit QED experiment. There, an ar- 




















7 


Fluence (cm'^) 


10 '^ 10 ” 10 ” 10 ” 10 ” 10 ” 10 ” 10 ” 


o 


10 “ 

10 “ 

10 " 

10 “ 

io“ 

10 ' 

10 ' 

10 “ 

10 " 

10 “ 

10 “ 

10 ' 



0.1 1 10 100 0.1 1 10 100 
Fluence (nm‘“) 


o 


Fluence (cm'“) 

10'“ lo” 10” lo” 



Luminescence intensity 
■ "S ->■ "I 

3/2 15/2 

• “H ->"l 

9/2 15/2 


Coupling strength 
□ m = 1/2 

'l6/2’ I 


(f) Model 1 (g) Model 2 







Figure 4. The luminescence intensity as a function of the implanted fluence is shown for the implanted at 600 K, (a) annealed and 
(b) non-annealed samples, also for implanted at 300 K, (c) annealed and (d) non-annealed samples. Fit curves are shown with 
the solid and dotted lines for three multiplets: '*F 7/2 —t'’li 5/2 (green), ^Hii /2 —>-'’li 5/2 (red) and '^ 83/2 —t“^ 115/2 (black). Solid lines 
relate to the fit with a single function, Eq. Dotted lines describe the fit with Model 1. (e) Application of RET model to the 
collective coupling strength, reported earlier by Probst et.al.^. (f) Schematic illustration of Model 1: after exceeding the critical 
fluence F 3 , ion-straggle increases while the concentration in the maximum (ml) remains the same, (g) Schematic illustration of 
Model 2: after the critical fluence F 3 is exceeded, ion-straggle increases as well as the concentration in the maximum (ml—l-ml*) of 
the ion distribution. 


ray of superconducting resonators is magnetically coupled 
to the room-temperature implanted sampl^i^, Fig.|^(e). 
The dependence of the collective coupling strength on the 
fluence is in agreement with Model 2. Accordingly, we 
obtain the modified fit-function for the collective coupling 
strength: 


9col-goJSresF ^ ( 6 ) 

V ^ 3 ^Fm dP- 

where gcoi is the collective coupling strength, go is the cou¬ 
pling strength per ion, Sres is the area under the resonator 
and Rpm is the Forster distance for the magnetic-dipole- 
magnetic-dipole interaction. Due to the lack of data-points 
in the region from 10^^ cm“^ to 10^^ cm“^, we cannot prop¬ 
erly fit the data by Model 1. Nevertheless, the solid line 
in Fig.|^(e) shows the closes correspondance to the data. 
Value of the Forster distance in this case is found to be 


equal to 2.76 nm, which is larger than 2.18 nm for the op¬ 
tical transitions. Thus, the magnetic dipole interaction oc¬ 
curs at larger distances between the ions. Nontheless, it 
explains the deviation of the coupling strength dependence 
on the number of ions at rather low fluenced^. 


V. CONCLUSION 

We studied the ion implantation process in Y 2 Si 05 crys¬ 
tals in combination with thermal annealing. We found that 
the proposed method produces stable and repetitive results 
for the ion activation in the crystals. The implantation 
yield for Erbium ions is estimated to be equal to 20%±5%. 
The luminescence studies suggest that the implanted ions 
experience more homogeneous crystal field as compared to 
the doped as-grown samples. However, the such homo- 

















































Figure 5. (a) Spectra of the samples annealed in different atmo¬ 
spheres. Spectrnm of each sample is normalized by the spectrum 
of doped as-grown sample, limp/las-grown- (b) An example of 
the relaxation-time measurement. 


geneous crystal field is accompanied by a larger number of 
defects, which are also stronger coupled to the Erbium ions. 
The luminescence intensity and collective coupling strength 
of implanted er ions reveals a non-linear dependence above 
the fluence of 10^‘^cm“^. Employing the resonance energy 
transfer theory, we developed two models describing the ob¬ 
served non-linear behaviours. The presented research paves 
the way towards a feasible technological process for fabri¬ 
cation of quantum-memory elements in hybrid circuits. 


Appendix A: Annealing procedures 


Several annealing procedures were applied for the post¬ 
implantation treatment. Rapid thermal annealing (RTA) 
was performed in N 2 atmosphere with a duration of 8 min¬ 
utes. Longer annealing lasted for 2 hours in two different 
atmospheres: air and Argon. Air was selected as oxygen- 
enriched atmosphere. Efficiency of the procedure is con¬ 
cluded from the similarity between optical spectra of the 
implanted and doped as-grown samples. Spectra of the 
implanted samples are normalized by the spectrum of the 
doped as-grown sample, as shown in Fig. [^(a). For the 
air-annealed samples (independently on the implantation 
temperature), normalized spectrum results in almost flat 
line. For the other samples, a lot of deviations are ob¬ 
served. Thus, the air-implanted sample has the local sym¬ 
metries most closely matching to the structure of the doped 
as-grown sample. 


Appendix B: Relaxation time 


Relaxation times as derives from pulsed measurements, 
when samples are exited with pulses of 50 /rs, see Fig.[^(b). 
Following decay of intensity is fit with an exponential func¬ 
tion, and thus the relaxation time values are derived. 


Appendix C: Resonance energy transfer. 


For the RET due to electric-dipole-electric-dipole inter¬ 
action, the quantum amplitude is given by: 

(Cl) 

where /r is the electric dipole operator and R is the inter¬ 
atomic distance, k is the wave vector, and A and B corre¬ 
spond to the donor and acceptor atoms. For the transitions 
at 560nm, the wave vector k = 1/560nm“^, and kR << 1 
at the interatomic distances below 20 nm, thus the k^R^ 
and k‘^R‘^ terms can be neglected. So, we have only the 
l/R® dependence: 

(C2) 


In a similar way, the quantum amplitude of the RET due 
to the magnetic-dipole-magnetic-dipole interaction is as fol¬ 
lows: 

,(C3) 

where m is the magnetic dipole operator. This similarly 
can be reduced to the dependence as above. For the 
electric-dipole-magnetic-dipole interaction, the main con¬ 
tribution comes from the term which is weaker than 

the an thus can be neglected in our calculations. 

The transfer rate function is derived from the depen¬ 
dence normalized by the Forster radius, which corresponds 
to the equal energy distribution between the radiative and 
the non-radiative energy transfer, thus the latter one equals 



(C4) 


and the final dependence for a single pair of atoms has the 
following form: 


Rad — R( 


1 -bF 


= Iu 


■ 


(C5) 


Appendix D: The fit function. 


The RET fit-function is used in the following form: 


7rl50^x 

^ 1-b 47rr|.x2/3dQ ’ 


(Dl) 


when the fluence x is smaller than the threshold fluence x < 
Fth- Above the threshold fluence, ion straggle d increases 
and also becomes a variable: 


7rl50^x 

^ 1 -b 47rr|.x2/3d2 


(D2) 
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In model 1, while the straggle increases, the concentration 
remains constant and the function reads: 


ttISO^x 


(D3) 
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7rl50^x 

^ ^°1+ 47rr|,x273d^„,^’ 
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